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Abstract 

We investigate the potential of 3 kiloTon-years(kTy) of KamLAND data to further con¬ 
strain the Am^ and tan^ 0 values compared to those presently allowed by existing KamLAND 
and global solar data. We study the extent, dependence and characteristics of this sensitivity 
in and around the two parts of the LMA region that are currently allowed. Our analysis with 
3 kTy simulated spectra shows that KamLAND spectrum data by itself can constrain Am^ 
with high precision. Combining the spectrum with global solar data further tightens the 
constraints on allowed values of tan^ 9 and Am^. We also study the effects of future neutral 
current data with a total error of 7% from the Sudbury Neutrino Observatory. We find that 
these future measurements offer the potential of considerable precision in determining the 
oscillation parameters (specially the mass parameter). 
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1 Introduction 


The Sudbury Neutrino Observatory (SNO) charged current (CC) and neutral current (NC) data on 
the measurement of solar neutrino flux has provided strong evidence for neutrino flavor conversion 
[1, 2] and it establishes the presence of flavor in the solar z/g flux at 5.5a level. The SNO and 
SK results together with the data from the radiochemical experiments Cl [3] and Ga [4] ^ single 
out the LMA solution based on MSW resonant matter conversion as the most probable solution 
of the solar neutrino puzzle [7]-[17]. In [9], for instance, this solution is characterized by best £t 
values of tan^ 6 = 0.41, Am^ = 6.06 x 10“^eV^ for the neutrino mixing parameters. Spectacular 
confirmation in favour of this solution, using terrestrial neutrino sources, has come recently from 
the KamLAND experiment in Japan [18]. The uniqueness of KamLAND [19, 20, 21] lies in its 
sensitivity to masses lying in the LMA region through its measurement of the reactor antineutrino 
energy spectrum. With 162 ton-years of data, KamLAND has already split the allowed LMA zone 
into two smaller sub-zones, the low-LMA with Am^ ~ 7.2 x 10“^ eV^ and the high-LMA zone 
with ~ 1.5 x 10"^ eV^ [22, 23]. 

In this paper, we consider a broader time-frame than spanned by the first results, and obtain 
projected results and conclusions that should be forthcoming from a 3 kTy exposure. With 
present fiducial volume this correspods to KL run of 8 - 10 years time. First we discuss the 
present constraints and next we discuss the future perspectives in the light of the present data. 
We determine the allowed areas in the Am^-tan^ 6 plane from the analysis of total KamLAND rate, 
KamLAND spectrum and examine the role of both of these in constraining the allowed regions 
by themselves as well as in conjunction with the global solar data. We simulate the projected 
KamLAND spectrum at several sample Am^ and tan^ 6 values, taken from the currently allowed 
regions, and attempt to delineate in detail the limits of sensitivity for KamLAND over a 3 kTy 
period. We study the dependence of the reconstituted parameter regions on the sample values of 
Am^ and tan^ 6 chosen and demonstrate that for spectrum simulated at points inside both the 
low-LMA and the high-LMA region the accuracy of reconstruction is quite high, leading to an 
excellent precision in determination of the oscillation parameters, especially Am^. 

Section 2 describes the salient features of the KamLAND detector, and the expected z7g reactor 
flux to which it is sensitive. In Section 3 we discuss the analysis procedure and results. Section 4 
summarises our conclusions. 


2 The KamLAND detector, the reactor flux of Electron 
Anti-Neutrinos and the Event Rate 

KamLAND [19] is a 1 kton liquid scintillator neutrino detector located at the earlier Kamiokande 
site in the Kamioka mine in Japan. It measures the z/g flux from 16 Japanese nuclear power 
reactors whose distances range from 80 km to 800 km. However 79% of the measured flux come 
from reactors situated at distances between 138 km to 214 km from the detector. The reaction 
that detects the z7g is the inverse beta decay z/g -|- p —> e"*" -|- n. The positrons are annihilated 

^For recent reviews on solar neutrino experiments see [5, 6]. 
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Table 1: Fitted Parameters a^, /c = 0,1, 2 for the reactor neutrino spectrum. The last row shows 
the energy released for each isotope per hssion. 


Isotope 


23yp^ 

238u 

^4iPu 

Oq 

0.870 

0.896 

0.976 

0.793 

Ql 

-0.160 

-0.239 

-0.162 

-0.080 

0-2 

-0.0910 

-0.0981 

-0.0790 

-0.1085 

(MeV) 

201.7 

205.0 

210.0 

212.4 


to produce two 7 rays. Neutron capture in the medium also generates a delayed 7 signal. The 
correlation between these two records an event grossly free from the backgrounds. 

The neutrino spectrum from the fission of a particular isotope j is conveniently parameterizable 
[21, 24] (in units of MeV“^ per fission) as 

dNi/dE^ = exp{ao + aiE^ + a2El). ( 1 ) 

Here j = 1, 2, 3,4, corresponding to the 4 isotopes 239p^ ^^^Pu which constitute the fuel. 

The htted values of the in the equation above are reproduced for completeness from [21, 24] 
in Table 1. 

In addition, each isotope has a characteristic energy released per fission, ej. These are also 
reproduced from [24] in Table 1. 

Table 2 (from [19]) gives the distances of the various reactors from the Kamioka mine which 
houses KamLAND, along with the maximum thermal power (in Giga-watts) of each i, the 

reactor index, which runs from 1 — 16. Also, we note that in principle, the power of each reactor 
varies over the year depending on demand, fuel composition, re-fuelling times etc. This dependence 
is averaged over for our purposes, and we assume that each reactor is in the running mode at 
maximum output 80% of the time. Using the above data, one may then write an expression for 
the spectrum from a given reactor i as, 

= E ^^^^k^^^dNl/dE, (2) 

j l^kJk^k 

Here /* is the fractional abundance of isotope j in reactor i at a given time. The time dependence 
is again averaged over, and for all the reactors we use for the abundance the values, 53.8% for 
32.8% for ^^®Pu, 7.8% for and 5.6% for ^^^Pu, as in [18]. Convenient units for Si are 
MeV“^sec“^, obtained by converting the power in Giga-watts into Mev per sec by multiplying by 
the factor 6.24 x 10^^. 

The other quantities needed to determine the event rate are the cross-section, the survival 
probability for anti-neutrinos and the number of free proton targets in the scintillator. The cross- 
section is given by 

Pe+Ee+. (3) 

mlfTn 
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Table 2: Reactor distances and power outpnts 


Reactor Site 

Distance di (km) 

Power ^^^^(Giga-watts) 

Kashiwazaki 

160 

24.6 

Ohi 

180 

13.7 

Takahama 

191 

10.2 

Hamaoka 

213 

10.6 

Tsuruga 

139 

4.5 

Shiga 

88 

1.6 

Mihama 

145 

4.9 

Fukushima-1 

344 

14.2 

Fukushima-2 

344 

13.2 

Tokai-II 

295 

3.3 

Shimane 

414 

3.8 

Ikata 

561 

6.0 

Genkai 

755 

6.7 

Onagawa 

430 

4.1 

Tomari 

784 

3.3 

Sendai 

824 

5.3 


Here / = 1.69 is the integrated Fermi fnnction for nentron /3-decay, me is the positron mass, 
Ee+ is the positron energy, Pe+ is the positron momentnm and = 886.7 secs is the nentron 
lifetime. The total visible energy {Eyis) corresponds to Ee+ + me, where Ee+ is the total energy 
of the positron and me the electron mass. The total positron energy is related to the incoming 
antineutrino energy E^, through the relation, Ee+ = E^ — {mn — mp) — Eyec, where mn — mp = 1.293 
MeV is the neutron-proton mass difference and E^ec is the average neutron recoil energy calculated 
here using [25]. 

The two-generation survival probability for the antineutrinos from each of the reactors i is 
given by 


T-./- -N . 9 /l.27Am^(iA ,,, 

Pi(z/e ^ z/g) = 1 - sm^26'sm^ f- — -j (4) 

where di is the distance of reactor i to KamLAND in km, E^ is in GeV and Am? is in eV^. Thus 
the total observed event-rate in KamLAND is given by N^l in sec“^, 

Nkl — V^p ^ J dEy{g J dEe+R(^Eyig, Ee+') x 

+ 1.293) (5) 
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where Np is the number of free protons in the hducial volume of the detector. KamLAND has 
declared 162 ton-year data corresponding to 408 tons hducial mass containing 3.46x10^^ free 
protons. The integrated total power P/AtiLF' is 254 Joule/cm^. The efficiency {rf) is 78.3%. 
R(E„i 5 ,Ee+) is the energy resolution function given by 


Ee+) 


1 f — Ef,+ 



with (To/Ee+ = 7.5%/^/E^ {E in MeV) [18, 26]. 

The ratio of the observed event-rate to the expected rate in KamLAND is dehned as 


Rkl — 




KL 


Kl 


( 6 ) 


where is the expected event-rate in KamLAND given by Eq.(5) with P{i'e ^ ^e) = L The 
hrst KamLAND results report a Rki = 0.611 ± 0.085 ± 0.041 [18]. 


3 Analysis and Results 


3.1 Current Data 


In Figure 1 we show the lines of constant KamLAND rate Rkl in the Am? — taffi 9 plane for the 
observed value of 0.611 and ± 3a limits. The best-ht point from global solar neutrino analysis and 
the corresponding 3a contour (dashed lines) in the LMA and LOW regions are also shown. The 
best-ht point [8, 9] of the global solar neutrino data predicts a KamLAND rate of 0.65 very close 
to the observed rate. We note that a particular observed Rkl may correspond to a wide range 
of KamLAND spectra. The reverse, of course, is not true since an observed KamLAND spectrum 
singles out a unique rate (within errors). 

We hrst do a statistical analysis with the KamLAND rate alone. For the rate we dehne the 


as 


Xkl 


/ T^expt Tytheory 
\^KL ~ ^KL 


2 


a 


2 


(7) 


where a = \Jalyg^ + ahati <^syst and astat being the total systematic and statistical error in the 
KamLAND data respectively. We take 6.42% systematic uncertainty along with the published 
statistical errors. To eliminate the geophysical background we take a visible energy threshold of 
2.6 MeV. [18, 26]. 

For a more complete statistical study, we next do a combined analysis of global solar data and 
the observed KamLAND rate. We use a combined function dehned as 


+ xIl ( 8 ) 

For the Xq (solar), we use the data on total rate from the Cl experiment, the combined rate 
from the Ga experiments (SAGE-I-GALLEX-I-GNO), the 1496 day data on the SK zenith angle 
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tan^ 0 


Figure 1: The isorate lines for the KamLAND detector corresponding to the observed rate and 
±3cr limits in the Am^ — tan^ 6 plane. Also shown by dotted line is the Sa allowed area from the 
global analysis of the solar neutrino data . The best-£t solution for the solar data is marked. 


energy spectrum and the combined SNO day-night spectrum. We dehne the function in the 
“covariance” approach as 

N 

xi=Y: («"■“ - - <“’) (9) 

where Ri are the solar data points, N is the number of data points (80 in our case) and 
is the inverse of the covariance matrix, containing the squares of the correlated and uncorrelated 
experimental and theoretical errors. For further details of our solar analysis we refer the reader 
to [9, 10]. 

In Figure 2 we draw the 90%, 95%, 99% and 99.73% C.L. allowed area in the LMA region 
from only KamLAND rate and combined solar-|-KamLAND rate analysis. Superimposed on that 
we show the 3(T allowed area from solar data alone. The current observed KamLAND rate is 
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Figure 2: The 90%, 95%, 99% and 99.73% C.L. contours in the LMA region. The left hand panel 
gives the allowed region from an analysis of only KamLAND rate. The right hand panel gives the 
allowed region from a combined solar+KamLAND rate analysis. The purple dashed line shows 
the 3cr only solar contour. 


in excellent agreement with the predicted KamLAND rate from the best-£t solution to the solar 
data. From the left-hand panel of Figure 2 we see that the allowed areas from KamLAND rate 
data are mostly consistent with the 3 ct allowed area from the global solar data, except from a 
small region around the bottom right-hand edge of the latter. This fact is again reflected in the 
global solar and KamLAND rate allowed areas shown in the right-hand panel. The inclusion of 
the KamLAND rate into the global analysis allows most of the regions allowed before, except for 
a small region around low Am^ and large tan^ 6 (the bottom right zone). 

Next, our aim is to see how far the allowed areas can be constrained with the inclusion of 
KamLAND spectrum data. The current KamLAND spectrum data are rather low on statistics 
and hence we consider a Poisson distribution for the spectral events. Thus the the current 
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Figure 3: The 90%, 95%, 99% and 99.73% C.L. contours in the LMA region. The left hand panel 
gives the allowed region from an analysis of only KamLAND spectrum data. The right hand panel 
gives the allowed region from a combined solar+KamLAND spectrum analysis. 


KamLAND spectrum is dehned as 


2 _ 
^KLspec / . 


Aj-exp 

2{XnNl^ - + 2N^^nn{ ^ 




+ 


- 1 )- 


a 


sys 


( 10 ) 


where agys is taken to be 6.42%, X^ is a normalization allowed to vary freely, and the sum is over 
the KamLAND spectral bins. 

The left hand panel in Figure 3 shows the allowed area from only KamLAND spectrum analysis. 
The right hand panel in Figure 3 gives the allowed region in the oscillation parameter space after 
including the global solar data. These contours are obtained by minimising 


— X© + XxLspec (11) 

Together the KamLAND and solar data are instrumental in narrowing down the parameter range 
by a large amount. At 99% C.L. the allowed LMA region is bifurcated into two parts - a low-LMA 
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lO 
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Figure 4: The la and Sa contours for the 3 kTy projected KamLAND spectrum analysis alone. 
The different panels are for the simulated spectrum at different hxed values of Am^ and tan^ 6 
shown by bold dots. For reference we also show by dashed lines the current 99% C.L. allowed 
areas from the global analysis of the solar and KamLAND data. 

zone around the best-£t point Am^ = 7.17 x 10“^ eV^ and tan^ 612 = 0.44 and a high-LMA zone 
around a second best-fit of Am^ = 1.49 x 10“^ eV^ and tan^ 612 = 0.43. At 3a the two regions 
merge. 


3.2 3 kton year data 

In this section we investigate the evolution of the allowed zones as KamLAND collects more 
statistics for its spectrum data. We do a projected analysis with 3 kTyr KamLAND spectrum 
simulated at few representative values of Am^ and tan^ 9. We randomize the generated spectra to 
take into account the possible fluctuations. We use these simulated spectra in a analysis and 
reconstruct the allowed regions in the Am^-tan^ 6 parameter space. Since statistics are expected 
to be large for the 3 kTy exposure, we consider a Gaussian distribution for the spectral events in 
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Figure 5: Ratio of oscillated to unoscillated events for the 3kTy simulated KamLAND spectrum 
for the different sets of Am^ and tan^0 corresponding to Figure 4. 


this case and dehne our function as, 

- A'f‘)(4)Kt(«T - X'”) (12) 

where {u‘fj)KL is the error correlation matrix, containing the statistical and systematic errors, where 
the systematic errors are assumed to be fully correlated among the energy bins. For the data we 
assume the same threshold of 2.6 MeV and the same energy binning as the present data. However 
the current KamLAND data quotes a conservative value of 6.42% for the systematic uncertainty. 
In the 3 kton year time span the systematic uncertainty is expected to reduce. The hducial volume 
uncertainty which contributes the most at present can go down after the collaboration installs a 
calibration arm. The cut systematics and energy reconstruction error can also reduce Keeping 

^We thank Prof. A.Suzuki, Prof. F. Suekane, Prof. S. Pakvasa and Prof. R. Svoboda for discussions on future 
systematic errors in KamLAND. 
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Figure 6: The la and 3cr contours for the combined analysis including the global solar data, the 
projected 3 kTy KamLAND data and the CHOOZ reactor data [27]. The different panels are 
for the simulated spectrum at values of Am^ and tan^ 9 of the corresponding panels of Figures 4 
and 5. The dashed line shows the presently allowed 3a contour from solar+KamLAND spectrum 
analysis. 


this in mind for 3 kTy data we use a optimistic assessment of 3% for the KamLAND systematic 
error, though we have checked that increasing the systematic uncertainty to 4% (which is the 
projected value quoted by the KamLAND collaboration) does not change the hnal conclusions. 

Figure 4 shows the reconstituted C.L. allowed contours (la and 3 a) in parameter space from 
a analysis using the simulated KamLAND spectrum alone. Marked by black dots are the 
points at which the spectra have been simulated. Each panel in this set of plots serves to identify 
regions in the presently allowed LMA region and demonstrates the constraining capabilities of 
the future KamLAND spectral data. Panels 1-5 are for spectrum simulated in the low-LMA zone 
while panels 6-9 are for spectrum simulated in the high-LMA zone. For the hrst three panels the 
spectrum is simulated at a Am^ corresponding to the low-LMA best-ht and three different 6s. 
A comparison of the allowed regions in these three panels show that tighter constraints on Am^ 
and tan^ 6 are associated with higher values of tan^ 9. For most of the cases where Am^ lies in 
the low-LMA region, the high-LMA can be ruled out at the 3a level. However for the panel 2 
which corresponds to the lowest 9 a small region still remains allowed at 3a. The panel 4 is for 
the spectrum simulated at the best-ht solar point, Am^ = 6.06 x 10“® eV^ and tan^ 6 = 0.41. As 
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Figure 7: Same as in Figure 6 but with a smaller future projected total uncertainty in the SNO 
NC rate of only 7%. 


we move our simulation point in the high-LMA region in panel 6, a higher Am^ region becomes 
allowed. No such region is obtained for the current best-£t point in the high-LMA region in panel 
7 as well as panel 8. Note that panel 8 is for the high-LMA best-£t Am^ but at an increased 6. 
Again in panel 9 we get some allowed regions at 3 ct at a high Am^ . We stress that since the 
relevant probability for KamLAND is the vacuum oscillation probability each of these panels will 
admit the mirror solution corresponding to 0 —> 7r/2 — 0, the so called dark zone, which we have 
not shown explicitly in this hgure. However we note that while very tight constraints are obtained 
for Am^ , the range of allowed values for tan^ 6 seems to be quite large in general for all the cases. 
In fact for all but panel 2, we see that maximal mixing is allowed atleast at the 3cr level, if not 
better. 

Figure 5 shows the distorsions in the simulated spectrum with ±lcr error bars at the Am^ and 
tan^ 6 corresponding to each of the panels of hgure 4. The spectra in panels 6 and 9 of Figure 5 
are simulated at a higher value of Am^ and hence this leads to faster oscillations as compared to 
the other panels, resulting in a degenerate solution appearing at a higher value of Am^ in hgure 

4. 

In hgure 6 we show the C.L. allowed regions from a combined analysis of the global solar and 
the 3 kTy KamLAND simulated spectrum data. The KamLAND spectra are simulated at the 
same points as in Figure 4 and 5. Also superimposed is the 99% C.L. combined solar -|- present 


12 


































KL contour. If we compare Figure 6 to Figure 4, we note that the contours get more constricted 
and the higher Am^ regions appearing in panels 2, 6 and 9 disappear. Inclusion of global solar 
data also restricts the allowed range in tan^ 6. The maximal mixing solution and the dark zone 
allowed by the KamLAND spectrum data get disfavoured with the inclusion of solar data in most 
of the panels, as the SNO data disfavours the maximal mixing solution. However for panels 3 
and 8 for which the simulation point is at a larger 6 the maximal mixing solution continues to be 
allowed at Sa, as these points correspond to spectra for which the KamLAND contribution to 
is much less at maximal mixing. 

In figure 7 we study the effect of reducing the SNO NC error in the analysis with 3 kTy 
KamLAND spectrum data. The current SNO NC data is due to neutron capture on deuteron the 
efficiency for which is 30%. SNO has added NaCl in the detector which has the effect of enhancing 
the neutron capture efficiency to 83%. With this the statistical error in SNO NC measurements is 
expected to reduce to about 5% [29]. The SNO collaboration are also improving their systematics. 
They expect to have less than 3% systematic uncertainty on neutron capture and would also 
reduce the uncertainties in their energy scale and energy resolution [29]. Keeping this in mind we 
assume an optimistic reduction in the current systematic error of 9% and we use a total SNO NC 
error of 7% [30] in place of the current 12% in figure 7. We do a y^-analysis of global solar data 
and the simulated KamLAND spectrum data. In plotting this figure, we use the total charged 
current and neutral current rates instead of the SNO spectrum data which implicitly assumes 
absence of spectral distortion for the resultant flux at SNO. Since we are in the LMA region 
we consider this to be a justified approximation above the SNO threshold of 5 MeV. This hgure is 
a most optimistic projection of how far the allowed ranges of and tan^ 6 can be constrained 
in the future. We hud that the reduction of the SNO NC error reduces the allowed range of 6 as 
compared to that in hgure 6. In particular in the panels 3,8 and 9 maximal mixing were allowed 
with 3 kTy KamLAND data but reducing the SNO NC error disfavours maximal mixing in these 
panels. 

In Table 4 we summarise the allowed ranges in Am^ and taffi 6 as obtained from analysis with 
solar + 3 kTy KamLAND spectrum data and solar (with sno NC 7%) and 3 kTy KamLAND spec¬ 
trum data. Note that for 3 kTy KamLAND spectrum we use an optimistic value of 3% for the 
KamLAND systematic error. We also give the spread [30] in each parameter where 

spread = ~ x 100 % (13) 

where V denotes the parameter Am? or taffi 9. 

We note that KamLAND has unprecedented sensitivity to Am^ and the 99% C.L. spread re¬ 
duces to 30% with the 0.162 kTy KamLAND data. With the inclusion of 3 kTy KamLAND data 
and improved systematics Am^ is determined with 6% precision. The 9 sensitivity on the other 
hand does not improve much with the inclusion of KamLAND data The neutrino energies cor¬ 
responding to the statistically signihcant region of the observed KamLAND spectra are close to 
the value required to make the probability (dehned in equation 4) one. Hence this region of the 
spectra has little sensitivity to theta and overall the sensitivity to the value of 9 is reduced. If we 
reduce the SNO NC error, the large values of the mixing angle are severely constrained and the 
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Data 

set 

used 

99% CL 
range of 

Am^ X 

lO-^eV^ 

99% CL 
spread 
of of 
Am? 

99% CL 
range 

in 

tan^ 6 

99% CL 
spread 
in 

tan^ 6 

only sol 

3.2 - 2.4 

76% 

.27- .75 

47% 

sol-|-162 Ty 

5.3 - 9.9 

30% 

.28 - .78 

47% 

sol-|-3 KTy 

6.8 - 7.7 

6 % 

.32 - .70 

37% 

sol(SNO NC 7%)+3 KTy 

6.9 - 7.7 

5.5% 

.33 - .60 

29% 


Table 3: The range of parameter values allowed. 

theta sensitivity improves. The probability relevant for SNO being the adiabatic MSW probability 
it goes as sin^ 9 giving it a greater sensitivity to constrain mixing angles 

4 Conclusions 

We have studied the capability of the future KamLAND data to constrain the mass and mixing 
parameters in the LMA region. We have simulated the KamLAND spectrum in the currently 
allowed LMA region for an exposure of 3 kTy. We hnd that with 3 kTyr exposure the Kam¬ 
LAND spectrum data is quite powerful in constraining Am^ giving unique allowed islands around 
the point at which the spectrum is simulated due to it’s high sensitivity to distortions driven by 
this parameter. We present a few cases where the spectra is somewhat flatter admitting multiple 
solutions. However this remaining ambiguity is removed with the inclusion of the solar data in 
the analysis irrespective of whether the spectrum is simulated at points belonging to low-LMA 
or high-LMA zone. Thus with 3 kton year exposure KamLAND can dehne Am^ very precisely 
to within 6%. The 6 sensitivity of KamLAND is not as good and after 3 kton year exposure 
and improved systematics the allowed spread in tan^ 6^ is ~ 37%. We hnd that for most of the 
simulated spectra the maximal mixing solution is disfavoured but may still get allowed at 3a if 
the true spectrum corresponds to higher value of tan^ 6. If however the SNO NC error is reduced 
then the maximal mixing gets disfavored in all the cases. 
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